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Abstract 

Correlations of like-sign pion pairs emerging from proton-antiproton collisions are 
analysed in the two-dimensional (qx, qx) decomposition of the three-momentum dif- 
ference q. While the data cannot be adequately represented by gaussian, exponen- 
tial, power-law or Edgeworth parametrisations, more elaborate ones such as Levy 
and an exponential with a cross term do better. A two-scale model using a hard cut 
to separate small and large scales may indicate a core that is more prolate than the 
halo. Consideration not only of the interference peak at small (qL,Qr), but also of 
the shape of the correlation distribution at intermediate momentum differences is 
crucial to understanding the data. 

PACS: 13.85.Hd, 13.87.Fh, 13.85.-t, 25.75.Gz 



1 Introduction 



The measurement of correlations between identical final-state particles, also 
commonly termed "Bose-Einstein correlations" or "HBT effect" after the fa- 
thers of intensity interferometry in astronomy, has become a valuable tool in 
the quest to understand the spacetime structure of high-energy collision pro- 
cesses [1,2]. The recently-coined term of femtoscopy highlights the possibilities 
of unraveling the interplay of kinematics and dynamics at the femtometer level 
based on the symmetrisation of identical bosons. 

The current concentrated effort at RHIC to quantify and understand ultrarel- 
ativistic nuclear collisions [3] relies extensively on comparisons with baseline 
scenarios constructed from the corresponding "trivial" hadron-hadron sample. 
In this context, the UAl experiment continues to be relevant and interesting, 



Preprint submitted to Elsevier Science 



7 February 2008 



even though data-taking at the CERN SPS has long ceased. Current experi- 
mental energies of 200 AGeV at RHIC are still below those available to UA1 
by a factor three, so that our results may also provide a window on possible 
energy dependencies of current investigations. 

One of the relevant theoretical frameworks is constructed in terms of a three- 
dimensional decomposition of the pair momentum difference, q = Pi — P2 
along the collision axis, the pair transverse momentum and the "side" direction 
orthogonal to both [4]. In this letter, we provide results on HBT analysis in 
terms of the simpler two-dimensional decomposition, defining in the usual way 
Ql — = |(cfz) z|, with z the beam direction, and qr = |q— qz,|- We analyse 
correlations between like-sign (LS) pion pairs in terms of the normalised second 
moment, 

MqL ' qT) -W(^)' (1) 

where p| lb (gL, qr) counts "sibling" LS pairs from the same event, while pf^qL, <?t) 
counts "reference" pairs made up through event mixing. 



2 Data sample, cuts and corrections 



Like-sign pion pairs from approximately 2.45 million minimum-bias events [5] 
measured by the UAl central detector in 1985 and 1987 were analysed. We 
applied the same single-track cuts used previously [6], including p±_ > 0.15 GeV 
\r]\ < 3 and 45° < \<f>\ < 135°. Good measurement quality and fitted track 
length £ > 30 cm were required. The sample contains mainly pions with an 
estimated 15% contamination of charged kaons and protons [7]. 

As in previous work [8,6], an angle cut was applied due to limited angular 
resolution in the central detector [9] , and pairs were required to have minimum 
four-momentum difference Q 2 = \J—(pi — P2) 2 > 0.0003 GeV 2 . 

Spurious "split-track" pairs strongly influence pair counts at small relative 
momenta. Besides eliminating these pairs with the same algorithm used in pre- 
vious analyses, we now also correct for the fact that this necessarily eliminates 
some physical like-sign track pairs. We determine correction factors for each 
(qi, qx) bin measured in the detector rest system and charged- multiplicity sub- 
sample (see below) by passing unlike-sign pairs through the same split-track 
algorithm. These correction factors are substantial at small (qL,qr), reach- 
ing 1.9 for low- multiplicity subsamples. Errors shown and quoted take into 
account the additional uncertainty introduced by the correction. 

We corrected for Coulomb repulsion by parametrising the Bowler Coulomb 
correction in the invariant momentum difference [10] with an exponentially 
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damped Gamov factor G(Q) [11], F coul (Q) = 1 + [G(Q) - 1] exp(-Q/Q eff ), 
with a best-fit value Q eS = 0.173 ± 0.001 GeV. 



Correlations from subsamples of fixed N, the observed charged multiplicity in 
the whole azimuth, were summed separately in the numerator and denomina- 
tor 

R (a) = Ejv^p 2 (q|A0 (2) 

2W J2 N PN(l-N^)p 1 ^p 1 (q\N)- {) 
The sum J2n Pn was implemented in terms of ten subsamples defined by multi- 
plicity N. The reference for each iV-subsample, pi<S>pi(q | N), was constructed 
by creating for each sibling event a set of fake events with the same multiplic- 
ity, using tracks randomly selected from the corresponding subsample track 
pool. The factor (1 — iV -1 ) = N(N— 1)/N 2 arises because the uncorrelated 
reference for fixed- N subsamples is not the poisson but the multinomial dis- 
tribution [12]. All HBT quantities are measured in the longitudinal co-moving 
system (LCMS) of the pion pair. 



3 Parametrisations 



We use the generic parametrisation R2{qL,qr) = 7 1 + A \S\ , where the 
multiplicative parameter 7 corrects for any remaining effect of the overall 
multiplicity distribution, A is the "chaoticity parameter" and S(qi,qT) corre- 
sponds to the Fourier-transformed source function. For \S\ 2 , we implemented 
the simple gaussian, the gaussian with cross-term [13], the simple exponential, 
and the exponential with cross-term parametrisations, 



\S\ 2 = exp(-Rlql-R 2 T qi), 

\S\ 2 = exp(-R 2 L q 2 L - R 2 T q 2 - 2R 2 LT q L q T ) , 

\S\ 2 = exp(-R L q L - R T qr) , 



\S\ = exp(-R L q L - R T q T - 2R LTy /q L q T ) , 
as well as a power-law parametrisation, 



Ri 
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1 + (R L q L y aL (Rrqr 



(3) 
(4) 
(5) 
(6) 



(7) 



Apart from the above parametrisations, nongaussian distributions can also be 
approached through Edgeworth expansions 1 [14,15], 



\S\ 2 = exp(-R 2 L q 2 L - R 2 T q 2 T ) J[ [l + K 4>d H A (V2 R d q d )/24 

d=L,T 



(8) 



1 Third-order cumulants must be zero due to the symmetry of R2(ql,Qt)- 
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or through Levy distributions [16], 



\s\ 2 = eM-Rhl-RW T T /2 , 

\S\ 2 = eM-Rhl - RWt - 2R 2 LT q L q T r /2 . 



(9) 
(10) 



All fits with the above parametrisations were performed over the whole region 
< q L , Qt < 0.5 GeV but omitting the bin (q L , q T ) < (0.02, 0.02) GeV which 
suffers from multiple large corrections and detector effects. Note that values for 
Rl and Rt cannot be compared between different parametrisations; note also 
that these parameters play the traditional role of "radius" of a distribution, 
and can be related to various source parameters [17], only for the gaussian 
cases (3)-(4). 



4 Analysis and Results 



4-1 The region of small (qL,Qr) 

Fig. 1 shows slices of i^C^Zo <?t) for fixed qx bins in the upper and fixed qi 
bins in the lower panels, together with the global fits based on Eqs. (3)- 
(10). The Gauss, Gauss- with-cross-term, and Edgeworth fits are practically 
indistinguishable, as Rlt in (4) is compatible with zero and cumulants k^l and 
k^t in (8) do not improve the fit. All three are equally bad, with x 2 /NDF ~ 
3.3. The simple exponential (5) (not shown) is even worse at 4.0. The cross- 
term exponential, however, has x 2 /NDF = 1.7. The power law (7) is a disaster 
at xVNDF = 16. 

Both Levy-based parametrisations (9)-(10) fare better in reproducing the 
strong peak observed in the data, with x 2 /NDF = 1.45 and 1.13 respectively. 
However, the fit parameters A, Rl, Rt, Rlt and a are strongly correlated so 
that there is no unique minimum and therefore no unique set of best-fit pa- 
rameter values. All that can be said with some confidence is that a is around 
0.20-0.23, that R L and R T are of order 10 3 fm with R L /R T = 1.6, and that 
R\ T is negative; we will return to this below. Omitting a second small- (g^, q?) 
bin from the Levy fits renders them even more unstable. This is hardly surpris- 
ing, as the above parameters collectively depend strongly on the exact shape 
of the peak in the very small (qL, qr) region, the very region that experimental 
measurement struggles to resolve. 
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4-2 Observations at intermediate (qL,qr) 



Fig. 1, revealing as it is, tells only part of the story as it focuses exclusively on 
the peak at small |q|. Intermediate scales, it turns out, are very important. 

In Fig. 2, the same data and some of the fits are shown as colour maps over 
the full set of bins, but with the peaks in the lower left corners truncated at 
R 2 = 1.9 in order to highlight structure at intermediate (qL,qr) scales. Panel 

(b) shows that the simple exponential (5) cannot possibly describe the data 
as its shape at intermediate scales is linear, while the data in (a) is elliptic. 
The Levy without cross term (9), Panel (e), matches the shape better. Panels 

(c) and (f), representing the exponential and Levy with cross terms included, 
most accurately reflect the data structure. The gaussian forms (4) and (8) have 
a shape similar to that in (e). From Panel (d), it is immediately apparent that 
the spectacular failure of the power law (7) stems from its hyperbolic shape. 
A "sum power law" R 2 — 7 [1 + (Rl qh) aL + [Rt <?t) _c * t ] fares even worse. 

Figure 2 also makes clear that the source — whatever its exact shape — is 
prolate, i.e. Rl > Rt must hold for all parametrisations. 



4-3 A two-scale approach. 



The strong peak seen in the data and the comparatively large values for x 2 
for the above parametrisations suggests that there may be two scales in the 
system. Borrowing the terms "core" and "halo" from the literature [18], one 
may try a simultaneous fit to a double-gaussian parametrisation [19,20], 

Rq(Ql, Qt) = 7 [l + Ac exp(-R 2 LC q 2 L - R 2 TC q T ) + \ H exp(-R 2 LH q 2 L - R TH q T ) 

(11) 

However, again there are too many parameters, so that there is no unique set 
of best-fit values. For this reason, we try a "two-scale" procedure, starting 
from the assumption that the two scales implicit in (11) can be separated 
by a hard cut in the data. First, we fit to the core gaussian only bins with 
momentum differences larger than a cutoff (qi > q cut or qr > q C ut), thereby 
fixing 7, Ac, Rlc an d Rtc- This core gaussian is then subtracted from all 
data, after which the remaining halo "data" with (qL,qr) < (<Zcut> ?cut) is fit 
with the halo gaussian 7 \ H exp(— R\ H q\ — RthQt) f° r fi xe d 7- Finally, the 
resulting core and halo fits are combined. 

This procedure must be tested for consistency. As shown in Fig. 3(a), the joint 
X 2 for both fits as a function of q cut = 0.02 nt,m is found to have a minimum in 
the interval g cut = 180-240 MeV. The joint best xV NDF = I- 28 for the two- 
scale model is comparable to those for the Levy fits and fit values are stable, so 



5 



that numbers can be quoted. The Rl, Rt, A and Rl/Rt values corresponding 
to the four smallest x 2 values are shown as filled points in Fig. 3(b)-(e). 
Averaging these numbers, we estimate Rlc — 0-75 ± 0.02 fm, R TC = 0.45 ± 
0.02 fm, R LH = 2.26 ±0.07 fm, and R TH = 1.89 ±0.10 fm, signalling a prolate 
core and a somewhat more spherical halo. We note that best-fit chaoticities 
A c = 0.23±0.01 and \ H = 0.92±0.04 are below the theoretical limit of 1, while 
the large intercept seen in the data itself (i?2(0, 0) > 3.7) and corresponding 
single-component chaoticities (A ~ 2.7) violate this limit. 

In Fig. 4, the upper lines show the combined two-scale best fit for q T and q^ 
slices, with the lower lines corresponding to the core. While there is substantial 
improvement over simple parametrisations especially at intermediate scales, 
data points in the peak remain consistently above even this fit. 

It must be emphasised that this two-scale procedure can work only if its 
assumptions are confirmed a posteriori. First, the resulting halo gaussian must, 
and does, become negligible at scales larger than q cnt for the ansatz to be valid. 
Also, the clear separation between the sizes of the core and halo radii seen in 
Fig. 3 do not contradict the assumption of the presence of two scales. For this 
data set, the two-scale model is consistent. 

4-4 Systematic errors and uncorrected data 

The above results have shown that the data appears to have a strong peak 
below 0.10 GeV which significantly exceeds all parametrisations tried. While 
it is tempting to conclude that the peak represents some physical effect, other 
possibilities must be checked. We hence conducted a survey of effects that 
various cuts and corrections have on R 2 . We find that the angle cut and Q 2 
cut and the restriction in azimuth have very little effect on the normalised 
moment and that any systematic error due to these is of the order of a few 
percent. 

However, the correction for the unwarranted removal of real LS track pairs by 
the split-track algorithm is large, and it is concentrated in precisely the region 
where the large peak occurs. Naturally, one must ask whether the entire excess 
of the data over various parametrisations reflects nothing but the correction 
itself. Taking the extreme approach of leaving the correction out altogether, 
we have repeated the entire analysis for the uncorrected data, which is shown 
in Fig. 5 together with the same set of parametrisations used in Fig. 1 plus 
the power law (7). 

It is immediately apparent that the uncorrected data peak cannot be de- 
scribed by gaussian parametrisations either (% 2 /NDF ~ 3.5). With regard to 
shape (not shown), the general elliptic and prolate form of the data remains 
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unchanged. The various parametrisations follow the pattern set by Fig. 2: 
again, the simple exponential and power law fail because they do not repro- 
duce the shape of the data; again, the Levy and exponential with cross term 
work better, with the latter faring best. While fit parameter values (e.g. of 
the Levy-fits with R L ~ 2.5-2.9 fm, R T ~ 1.6-1.8 fin and a ~ 0.7) differ from 
those of the corrected data, the conclusions regarding the success of particular 
parametrisations are hence independent of the split-track correction. 



5 Discussion and Conclusions 

The present results on HBT correlations extend UAl observations to multi- 
dimensional correlations for the first time, which, at 630 GeV, represent the 
highest CMS energy at which this has been done. Two issues stand out: the 
strength of the peak at small momentum differences, and the importance of 
shape at intermediate scales. 

With a strong peak at small qL,£[T < 0.1 GeV, our high-statistic data, with 
and without corrections, rules out parametrisations based on single gaussians 
and their derivatives. While some other parametrisations are strongly peaked, 
none of those tried, not even the Levy and two-scale cases, reproduces the peak 
data convincingly. This may hint that "something else is going on" , be it the 
influence of jets, resonances, clustering effects or some other unknown factor. 
The two-scale results may hint that the "halo" could be due to short-lived 
resonances such as the p. 

We stress that it is unlikely that the strong peak seen in our data at small 
momentum differences is due to bias. First, the peak persists even without 
the Coulomb or split-track corrections. Second, the nongaussian behaviour 
seen in this paper is in line with UAl one-dimensional correlation structures 
seen earlier in the form of nonzero higher-order cumulants [7,21] and in the 
power law in Q 2 [8,22], even though the data had not been corrected in the 
way it has been here. Third, the presence of unidentified kaons and protons 
in the sample imply that the real peak should exceed the one shown here, so 
that the present results are conservative. It should be noted that a number 
of other hadronic experiments [23,24,25] have also seen significant deviations 
from gaussian behaviour at small |q|. 

Besides the structure of the peak seen at small scales, the shape of the distribu- 
tion at intermediate scales is seen to provide valuable additional information. 
Fig. 2 shows that the data has without doubt an elliptic and prolate shape, 
the latter in the sense that the distribution in momentum space is narrower 
in the longitudinal than in the transverse direction. This is confirmed also 
within the two-scale method and for the uncorrected data. Significantly, these 
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conclusions are independent of the various parametrisations and corrections. 

Second, plots of shape (whether in colour or as contour lines) help to constrain 
possible parametrisations: whereas in the one- dimensional analysis of UAl 
data [8], the exponential and power-law parametrisations were found to be 
superior to the gaussian, their extension to two dimensions fails badly because 
their contour lines are straight lines or even hyperbolic. 

Considering shape in terms of different parametrisations is complementary 
to decompositions into cartesian and spherical harmonics [26,27], which can 
be expected to work best for near-gaussian data. We note, however, that the 
interrelationship between shape and algebraic form of a parametrisation is 
less than obvious: elliptic shape is exhibited not only by bilinear forms such 
as Eqs. (3)-(4) and (9)-(10) but also, surprisingly, by the exponential-with- 
cross-term, whose success depends strongly on the sign of the cross term. For 
both exponential (6) and Levy with cross-term (10), our data clearly prefers 
negative values for the cross term (Rlt and R\ T respectively); positive values 
result in more hyperbolic shapes, in conflict with the data. 

Their fit instability nonwithstanding, the Levy parametrisations, which have 
the right shape and an acceptable peak, appear promising. Taking these at 
face value would imply that the corresponding source falls off with a power-law 
tail at large distances [16,28,29]. 
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Fig. 1. Upper panels: R2(qL, Or) data and best fits, shown left to right for slices with 
fixed q T = 0.00-0.02, 0.02-0.04, 0.04-0.06 and 0.06-0.08 GeV bins. Solid lines: Levy 
fit (10); dashed: exponential with cross term; dash-dotted: Gauss/Edgeworth. Lower 
panels: R2(ql,Qt) and the same fits for corresponding fixed-gx slices. Throughout 
this letter, fits are performed over all 624 data points, not just the data points, slices 
and intervals shown in the figures. 
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Fig. 2. Shapes of R2(qL,QT) data and fits at intermediate scales. Panels show (a) 
UAl data, (b) exponential, (c) exponential with cross term (d) power-law, (e) Levy, 
(f) Levy with cross term. Plots are truncated vertically to R2 < 1.9 to highlight 
structure at intermediate scales. 
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Fig. 3. Dependence on q cut - (a) Combined \ 2 f° r core an d halo fits as a 
function of g cu t = (0.02 n^) GeV. The four lowest \ 2 values correspond to 
q cut = 0.18-0.24 GeV and x 2 /NDF = 1.28. (b)-(e): Dependence of parameter val- 
ues on q cut . Filled points correspond to the four smallest \ 2 m ( a )- 
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Fig. 4. R2(qL,QT) and fit to core (lower lines) and combined core-halo (upper lines), 
shown in the upper panels for slices with fixed qi = 0.00-0.02, 0.02-0.04, 0.04-0.06 
and 0.06-0.08 GeV and in the lower panels for corresponding fixed-g^ slices. 
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Fig. 5. Same as Fig. 1 but for uncorrected data and including power law (dotted 
line). Note the different vertical scale. 
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